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Abstract 

The solar neutrino problem has persisted for almost three decades. Recent results from 
Kamiokande, SAGE, and GALLEX indicate a pattern of neutrino fluxes that is very 
difficult to reconcile with plausible variations in standard solar models. This situation is 
reviewed and suggested particle physics solutions are discussed. A summary is given of the 
important physics expected from SNO, SuperKamiokande, and other future experiments. 
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1. INTRODUCTION 

Thirty years ago, in the summer of 1965, the Homestake Mining Company completed 
the excavation of the 30 x 60 x 32 ft cavity that was to house the 100,000-gallon chlorine 
detector proposed by Ray Davis Jr. and his Brookhaven National Laboratory collabora- 
tors. Three years later Davis, Harmer, and Hoffman (1968) announced the results of their 
first two detector runs, an upper bound on the solar neutrino flux of 3 SNU (1 SNU = 
10~^^ captures/ target atom/ sec). The accompanying theoretical paper by Bahcall, Bah- 
call, and Shaviv (1968) found a rate of 7.5 ± 3.3 SNU for the standard solar model. This 
discrepancy persists today, augmented by nearly three decades of data from the Home- 
stake experiment and by the new results from the Kamiokande and GALLEX/SAGE 
detectors. The purpose of this review is to summarize the current status of the solar 
neutrino problem and its possible implications for physics and astrophysics. 

A remarkable aspect of the solar neutrino problem is that it has both endured and 
deepened. The prospect of quantitatively testing the theory of main-sequence stellar 
evolution provided much of the original motivation for measuring solar neutrinos: solar 
neutrinos carry, in their energy distribution and flux, a precise record of the thermonu- 
clear reactions occurring in the sun's core. Our understanding of the atomic and nuclear 
microphysics governing stellar evolution - nuclear reaction rates, radiative opacities, and 
the equation of state - has progressed significantly since the 1960s. The development 
of helioseismology has provided a new tool for probing the solar interior. Finally, we 
better understand our sun in the context of other stars, the observations of which have 
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helped to define the envelope of possibilities for diffusion, mass loss, magnetic fields, etc. 
This progress has tended to increase our confidence in the standard solar model. At the 
same time, with the new results from Kamiokande and the gallium detectors, a pattern 
of neutrino fluxes has emerged that is more difficult to reconcile with possible variations 
in that model. 

The solar neutrino problem has deepened because of the discovery of the Mikheyev- 
Smirnov-Wolfenstein (MSW) mechanism: the sun has the potential to greatly enhance 
the effects of neutrino mixing. If neutrino oscillations prove to be the solution to the 
solar neutrino problem, this will force modifications in the standard model of electroweak 
interactions, which accommodates neither massive neutrinos nor their mixing. This new 
physics would have implications for a variety of problems in astrophysics, including the 
missing mass puzzle and the formation of large-scale structure. 

There are several sources that the interested reader can consult for additional in- 
formation. The most comprehensive treatment is that given by Bahcall in his book 
"Neutrino Astrophysics" (1989). The appendix of this book reprints a delightful histor- 
ical perspective of the development of the solar neutrino problem (Bahcall and Davis 
1982). The review by Bowles and Gavrin (1993) provides an excellent discussion of 
the Homestake, Kamiokande, and SAGE/GALLEX experiments, as well as detectors 
now under construction or development. The proceedings of the Seattle Solar Model- 
ing Workshop (Balantekin and Bahcall 1995) and the Homestake Conference (Cherry, 
Lande, and Fowler 1985) contain many studies of the standard and nonstandard solar 
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models and of the nuclear and atomic physics on which such models depend. The theory 

of matter-enhanced neutrino oscillations has been reviewed recently by Mikheyev and 

Smirnov (1989). 

2. SOLAR MODELS 

2.1 The Standard Solar Model 

Solar models trace the evolution of the sun over the past 4.6 biUion years of main 
sequence burning, thereby predicting the present-day temperature and composition pro- 
files of the solar core that govern neutrino production. Standard solar models share four 
basic assumptions: 

* The sun evolves in hydrostatic equilibrium, maintaining a local balance between the 
gravitational force and the pressure gradient. To describe this condition in detail, one 
must specify the equation of state as a function of temperature, density, and composition. 

* Energy is transported by radiation and convection. While the solar envelope is con- 
vective, radiative transport dominates in the core region where thermonuclear reactions 
take place. The opacity depends sensitively on the solar composition, particularly the 
abundances of heavier elements. 

* Thermonuclear reaction chains generate solar energy. The standard model predicts 
that over 98% of this energy is produced from the pp chain conversion of four protons 
into ^He (see Figure 1) 

4p He + 2e+ + 2i/e (1) 
with the CNO cycle contributing the remaining 2%. The sun is a large but slow reac- 
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tor: the core temperature, Tc ~ 1.5 • 10^ K, results in typical center-of-mass energies for 
reacting particles of ~ 10 keV, much less than the Coulomb barriers inhibiting charged 
particle nuclear reactions. Thus reaction cross sections are small, and one must go to sig- 
nificantly higher energies before laboratory measurements are feasible. These laboratory 
data must then be extrapolated to the solar energies of interest. 

* The model is constrained to produce today's solar radius, mass, and luminosity. 
An important assumption of the standard model is that the sun was highly convective, 
and therefore uniform in composition, when it first entered the main sequence. It is 
furthermore assumed that the surface abundances of metals (nuclei with A > 5) were 
undisturbed by the subsequent evolution, and thus provide a record of the initial solar 
metallicity. The remaining parameter is the initial ^He/H ratio, which is adjusted until 
the model reproduces the present solar luminosity after 4.6 billion years of evolution. 
The resulting ^He/H mass fraction ratio is typically 0.27 ± 0.01, which can be compared 
to the big-bang value of 0.23 ± 0.01 (Walker et al. 1991). Note that the sun was formed 
from previously processed material. 

The model that emerges is an evolving sun. As the core's chemical composition 
changes, the opacity and core temperature rise, producing a 44% luminosity increase since 
the onset of the main sequence. The neutrino flux, the most temperature-dependent 
component, proves to be of relatively recent origin: the predicted flux increases exponen- 
tially with a doubling period of about 0.9 billion years. The equilibrium abundance and 
equilibration time for ^He are both sharply increasing functions of the distance from the 
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solar center. Thus a steep ^He density gradient is established over time. 

The principal neutrino-producing reactions of the pp chain and CNO cycle are sum- 
marized in Table 1. The first six reactions produce 13 decay neutrino spectra having 
allowed shapes with endpoints given by E™^ . Deviations from an allowed spectrum 
occur for ®B neutrinos because the ®Be final state is a broad resonance (Bahcall and 
Holstein 1986); much smaller deviations occur because of second-forbidden contributions 
to the decay. The last two reactions produce line sources of electron capture neutrinos, 
with widths ~ 2 keV characteristic of the solar core temperature (Bahcall 1993). The 
resulting solar neutrino spectrum is shown in Figure 2. 

Measurements of the pp, ^Be, and ^B neutrino fiuxes will determine the relative con- 
tributions of the ppl, ppll, and ppIII cycles to solar energy generation. As the discusion 
below will illustrate, this competition is governed in large classes of solar models by a 
single parameter, the central temperature Tg. The fiux predictions of two standard mod- 
els, those of Bahcall and Pinsonneault (1992) and of Turck-Chieze and Lopez (1993), are 
included in Table 1. 

2.2 Uncertainties in Standard Solar Model Parameters 

Careful analyses of the experiments that will be described in Section 3 indicate that 
the observed solar neutrino fiuxes differ substantially from standard solar model (SSM) 
expectations (White, Krauss, and Gates 1993; Parke 1995; Hata and Langacker 1994): 

0(pp) ~ 0.9(/)^^^(pj)) 
0(^Be) ~ 
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O_43 0SSM(8b)_ 



(2) 



Reduced Be and neutrino fluxes can be produced by lowering the central temperature 
of the sun somewhat. However, such adjustments, either by varying the parameters of 
the SSM or by adopting some nonstandard physics, tend to push the 0(^Be)/0(^B) ratio 
to higher values rather than the low one of Eq. (2), 



Thus the observations seem difficult to reconcile with plausible solar model variations. 

It is apparent that the rigor of this argument is the crucial issue: how quantitative is 
the tracking of fluxes and flux ratios with Tc, what variations can exist in models that 
produce the same T^. but differ in other respects, and how significant are the results in 
Eq. (2) when the statistical and systematic errors of the experiments arc taken into 
account? These questions have motivated a number of careful examinations (and clever 
presentations) of solar model uncertainties. 

SSM uncertainties include the reaction cross sections for the pp chain and CNO cycle, 
the opacities, the deduction of heavy element abundances from solar surface abundances, 
the solar age and present day luminosity, and the equation of state. Modelers occasionally 
adopt different "best values" and associated errors for these parameters. Bahcall and 
Pinsonneault (1992) have argued that the differences among solar models are almost 
always attributable to parameter choices, and not to disagreements about the underlying 
physics. 

While a detailed summary of standard model uncertainties would take us well beyond 
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the limits of this review, a quahtative discussion of pp chain uncertainties is appropriate. 
This nuclear microphysics has been the focus of a great deal of experimental work, as 
well as the source of some contention among modelers. The pp chain involves a series of 
nonresonant charged-particle reactions occurring at center-of-mass energies that are well 
below the height of the inhibiting Coulomb barriers. As the resulting small cross sections 
preclude laboratory measurements at the relevant energies, one must extrapolate higher 
energy measurements to threshold to obtain solar cross sections. This extrapolation is 
often discussed in terms of the astrophysical S-factor (Fowler 1984; Burbidge, Burbidge, 
Fowler, and Hoyle 1957), 



ity of the coUiding particles. This parameterization removes the gross Coulomb effects 
associated with the s-wave interactions of charged, point-like particles. The remaining 
energy dependence of S(E) is gentle and can be expressed as a low-order polynomial 
in E. Usually the variation of S(E) with E is taken from a direct reaction model and 
then used to extrapolate higher energy measurements to threshold. The model accounts 
for finite nuclear size effects, strong interaction effects, contributions from other partial 
waves, etc. As laboratory measurements are made with atomic nuclei while conditions 
in the solar core guarantee the complete ionization of light nuclei, additional corrections 
must be made to account for the different electronic screening environments. 

Among the pp chain reactions, the one presently most controversial is ^Be(p,7)®B. 




(4) 



where rj — ^f^" , with a the fine structure constant and P 



v/c the relative veloc 
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The two experiments performed at the lowest energies (Kavanagh et al. 1969; Fihppone 
et al. 1983) find very similar energy dependences but disagree in overall normalization 
by about 25%. A similar normalization disagreement exists among higher energy data 
sets (see Langanke 1995; Parker 1968; Kavanagh et al. 1969; Vaughn et al. 1970). 

Theory predicts a threshold energy dependence for S(E) that rises gently as E — > 0; 
at sufficiently low energies, when the capture occurs at nuclear separations well outside 
the range of nuclear forces, the behavior is fixed (Christy and Duck 1961; Williams 
and Koonin 1981). The data are consistent with theory, though their accuracy at low 
energies (120-200 keV) is insufficient to independently determine the rise at low E (see 
Figure 3). When the theoretical energy dependence (Johnson et al. 1992) of S(E) is used 
to extrapolate the Kavanagh et al. and Filippone et al. data, one finds S(0) ~ 25 eV 
barns and ~ 20 eV barns, respectively. While most calculations have predicted values 
for S(0) closer to the Kavanagh value, lower values have also been obtained (see, e.g., Xu 
et al. 1994) and general arguments have been given that the theoretical normalization 
uncertainty could be as much as a factor of two (Barker 1980). Thus it does not appear 
that theory can distinguish between competing experimental normalizations for S(0). In 
the most recent reexamination of the data, Johnson et al. (1992) determined 817(0) = 
22.4 ± 2.1 eV barns, from a weighted average of six data sets. 

Very recently Motobayashi et al. (1994), by measuring the breakup of in the 
Coulomb field of ^°^Pb, deduced a preliminary value of 817(0) = 16.7 ± 3.2 eV barns, 
which would favor the Fihppone result. It has been argued that the extraction of 817(0) 
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from the Coulomb breakup cross section is complicated by the three-body Coulomb 
effects in the outgoing channel and by larger E2 contributions (Langanke and Shoppa 
1994), though this claim has been disputed (Gai and Bertulani 1994). (Transverse electric 
multipoles of the electromagnetic current operator of rank J are denoted EJ.) 

The SSM calculations of Bahcall and Ulrich (1988) and Turck-Chieze et al. (1988) 
adopted 817(0) values of 24.3 ± 1.8 eV b and 21 ± 3 eV b, respectively. This was the 
most important contribution to the differences in their flux predictions. Their recently 
updated calculations (Bahcall and Pinsonneault 1992; Turck-Chieze and Lopez 1993) 
both adopt the Johnson et al. central value, though Turck-Chieze and Lopez assign a 
larger error bar to this result because of the concern that the weighted average of Johnson 
et al. does not fully reflect the systematic disagreement illustrated in Figure 3. 

While there has been some movement to lower values of 817(0), neither the best 
value nor appropriate error is likely to be decided without additional measurements. The 
resulting reduction in 0(^B) because of the weaker branch to the ppIII cycle reduces the 
discrepancy between the 88M predictions and the ^^Cl results. This reduction does not 
resolve the solar neutrino problem because it also leads to an increase in the predicted 
ratio 0(''Be)/0(^B), exacerbating the most puzzling aspect of Eq. (2). 

8mall differences between the Bahcall and Pinsonneault (BP) and Turck-Chieze and 
Lopez (TCL) calculations also exist in the low-energy extrapolation of the ^He-^He data 
of Krauss et al. (1987). BP adopt the value of Parker and Rolfs (1991), who correct for 
small electron screening contributions that are believed to enhance the laboratory cross 
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sections obtained at the lowest energies. The result is 833(0) = 5.0 ± 0.3 MeV b. TCL 
adopted the central value suggested by Assenbaum et al. (1987) of 5.24 MeV b. The 
earher Bahcall and Ulrich (1988) and Turck-Chieze et al. calculations had used central 
values of 5.15 and 5.57 MeV b, respectively, the latter from a direct extrapolation of the 
Krauss et al. data. The somewhat lower values adopted in the most recent calculations 
strengthen the branches to the ppll and ppIII cycles. 
2.3 Standard Solar Model Neutrino Fluxes 

There are some additional minor differences in the BP and TCL SSMs, e.g., in the 
solar lifetime, in 834(0), in the treatment of plasma effects on Thompson scattering, and 
in the composition. In addition, the BP calculation differs from TCL, and from the 
earher Bahcall and Ulrich (BU) SSM, by the inclusion of helium diffusion, which now 
is the largest contributor to the differences in the resulting flux predictions. The fluxes 
for the TCL and BP (with and without hehum diffusion) SSMs are given in Table 1. It 
is reassuring that equivalent calculations (TCL and BP without helium diffusion) differ 
by less than 10% in their predictions of the temperature-dependent 4>{^B). This flux is 
increased by 12% when helium diffusion is included. 

In the remainder of this review, citations to the BP SSM will refer to their model 
with helium diffusion. 

More important than the "best values" of the fluxes are the ranges that can be 
achieved by varying the parameters of the SSM within plausible bounds. In order to 
take into account the correlations among the fluxes when input parameters are varied. 
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BU constructed 1000 SSMs by randomly varying five input parameters, the primordial 
heavy-element-to-hydrogen ratio Z/X and S(0) for the p-p, ^He-^He, ^He-^He, and p- 
^Be reactions, assuming for each parameter a normal distribution with the mean and 
standard deviation used in their 1988 study. (These were the parameters assigned the 
largest uncertainties.) Smaller uncertainties from radiative opacities, the solar luminosity, 
and the solar age were folded into the results of the model calculations perturbatively, 
using the partial derivatives of the BU SSM (Bahcall and Haxton 1989). 

The resulting pattern of ^Be and ^B flux predictions is shown in Figure 4. The 
elongated error eUipses indicate that the fluxes are strongly correlated. Those variations 
producing 0(^B) below O.80^^^^(®B) tend to produce a reduced 0(^Be), but the reduction 
is always less than 0.8. Thus a greatly reduced 0(^Be) cannot be achieved within the 
uncertainties assigned to parameters in the SSM. 

A similar exploration, but including parameter variations very far from their preferred 
values, was carried out by Castellani, Degl'Innocenti, Fiorentini, Ricci, and collaborators 
(1994 and 1995), who displayed their results as a function of the resulting core tem- 
perature Tp. The pattern that emerges is striking (see Figure 5): parameter variations 
producing the same value of Tc produce remarkably similar fluxes. Thus provides 
an excellent one-parameter description of standard model perturbations. Figure 5 also 
illustrates the difficulty of producing a low ratio of 0(^Be) / 0(^B) when Tc is reduced. 

The BU 1000-solar-model variations were made under the constraint of reproducing 
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the solar luminosity. Those variations show a similar strong correlation with Tc 

cj>{pp)^T;'-' 0rBe)cxT« 0f B) oc Tf . (5) 

Figures 4 and 5 are a compelling argument that reasonable variations in the parameters 
of the SSM, or nonstandard changes in quantities like the metallicity, opacities, or solar 
age, cannot produce the pattern of fluxes deduced from experiment (Eq. (2)). This 
would seem to limit possible solutions to errors either in the underlying physics of the 
SSM or in our understanding of neutrino properties. 
2.4 Nonstandard Solar Models 

Nonstandard solar models include both variations of SSM parameters far outside the 
ranges that are generally believed to be reasonable (some examples of which are given in 
Figure 5), and changes in the underlying physics of the model. The solar neutrino problem 
has been a major stimulus to models of the second sort. But there are other observations 
that suggest nonstandard processes could take place in the sun. The depletion of ''Li in 
our sun (by a factor of 200) and in other solar-like stars is not understood. In a recent 
summary, Charbonnel (1995) argued that the depletion occurs on the main sequence 
after ~ 10^ years and increases with time. The SSM predicts no solar lithium depletion: 
the base of the convective zone is too shallow to reach temperatures where lithium can 
be burned. While lithium depletion in lower mass stars (~ O.9M0) can be produced 
by lowering the low-temperature opacities, the depletion in heavier stars would seem to 
argue for diffusion or mixing mechanisms not yet incorporated into the SSM (Charbonnel 
1995). 
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It is far from clear that the hthium depletion problem will require changes in the 
SSM that will affect the deep, neutrino-producing regions of the sun. On the contrary, 
Charbonnel argues that the observed ^^C/^^C ratios in low- mass red giants requires a 
diffusion coefficient that decreases rapidly with depth. 

Many nonstandard models were constructed to produce a reduction in Tc of about 
5%, as this would have accounted for the low counting rate found in the Homestake 
experiment. The suggestions included models with low heavy element abundances ("low 
Z" models), in which one abandons the SSM assumption that the initial heavy element 
abundances are those we measure today at the sun's surface; periodically mixed solar 
cores; models where hydrogen is continually mixed into the core by turbulent diffusion 
(Schatzman and Maeder 1981) or by convective mixing (Ezer and Cameron 1988; Shaviv 
and Salpeter 1968); and models where the solar core is partially supported by a strong 
central magnetic field (Abraham and Iben 1971; Bartenwerfer 1973; Parker 1974; Ulrich 
1974) or by its rapid rotation (Demarque, Mengel, and Sweigart 1973), thereby relaxing 
the SSM assumption that hydrostatic equilibrium is achieved only through the gas pres- 
sure gradient. A larger hst is given by Bahcall and Davis (1982). To illustrate the kinds 
of consequences such models have, two of these suggestions are discussed in more detail 
below. 

In low-Z models (Iben 1969; Bahcall and Ulrich 1971) one postulates a reduction in 
the core metallicity from Z ~ 0.02 to Z ~ 0.002. This lowers the core opacity (primarily 
because metals are very important to free-bound electron transitions), thus reducing Tc 
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and weakening the ppll and ppIII cycles. The attractiveness of low-Z models is due in part 
to the existence of mechanisms for adding heavier elements to the sun's surface. These 
include the infall of comets and other debris, as well as the accumulation of dust as the 
sun passes through interstellar clouds. However, the increased radiative energy transport 
in low-Z models leads to a thin convective envelope, in contradiction to interpretations 
of the 5-minute solar surface oscillations (Rhodes, Ulrich, and Simon 1977; Christensen- 
Dalsgaard and Gough 1980). A low He mass fraction also results. Furthermore, Michaud 
(1977), noting that diffusion of material from a thin convective envelope into the interior 
would deplete heavy elements at the surface, questioned whether present abundances 
could have accumulated in low-Z models. Finally, the general consistency of solar heavy 
element abundances with those observed in other main sequence stars makes the model 
appear contrived. 

Models in which the solar core (~ 0.2 Mq) is intermittently mixed (Dilke and Gough 
1972; Fowler 1972; Opik 1953) break the standard model assumption of a steady-state 
sun: for a period of several million years (approximately the Kelvin-Helmholtz time for 
the core) following mixing, the usual relationship between the observed surface lumi- 
nosity and rate of energy (and neutrino) production is altered as the sun burns out of 
equilibrium. Calculations show that both the luminosity and the neutrino flux are 
suppressed while the sun relaxes back to the steady state. 

Such models have been considered seriously because of instabilities associated with 
large gradients in the ^He abundance. In the higher temperature central regions of the 
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sun the pp chain reaches equihbrium quickly. The dominant ppl cycle production of ^He 
is controlled by the p+p reaction rate which varies as X^T^, where X and T are the local 
H mass fraction and temperature, respectively, while it is destroyed by the ^He+^He 
reaction at a rate proportional to X3 T^^. Thus, the core ^He abundance X3 oc and 
rises steeply with increasing radius until a point where '^He equilibrium has not yet been 
reached in the pp chain. This is approximately the peak of the ^He abundance (at r ~ 
O.3r0 is the SSM): beyond this point the abundance declines quickly. (See Bahcall 1989, 
Figure 4.2.) This profile is unstable under finite amphtude displacements of a volume to 
smaller r: the energy released by the increased ^He burning at higher T can exceed the 
energy in the perturbation. 

Dilke and Gough (1972), in the "solar spoon," proposed that the increased ^He burn- 
ing would also produce a linear instability of low-order, low-degree gravity modes which 
they postulated could trigger large-scale mixing of the solar core. Under small oscil- 
lations, the enhanced burning in a volume element displaced downward increases the 
bouyant restoring force, leading to greater ascending velocities than descending ones. 
The question is whether this instability survives damping mechanisms such as radiative 
difi^usion, turbulent convection, and couplings to higher order g-modes. As Merryfield 
(1995) recently summarized, the possibility of an instability is still open, though no calcu- 
lation has included all of the damping mechanisms thought to be important. He is far less 
optimistic about the possibility that the instability would drive large-scale, intermittent 
mixing (Merryfield, Toomre, and Gough 1990). Theory suggests that the amphtude of 
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such oscillations would saturate at velocities of ~ 10 km/s; but observation seems to rule 
out such large amplitude, low-order, low-degree g-modes. The expected surface velocities 
of such modes would be comparable to the core velocities, exceeding the helioseismology 
bound of ~ 10 cm/s substantially. 

These discussions of two of the more seriously explored nonstandard models illustrate 
how changes motivated by the solar neutrino problem often produce other, unwanted 
consequences. Both examples underscore the growing importance of helioseismology as 
a test of the SSM and as a constraint on its possible variations. 

Figure 6 is an illustration by Hata (1995) of the flux predictions of several nonstandard 
models, including a low-Z model consistent with the ^''Cl results. As in the Castellani et 
al. exploration, the results cluster along a track that defines the naive Tc dependence of 
the 0(^Be)/0(^B) ratio, well separated from the experimental contours. 

Of course, one cannot rule out a nonstandard solar model solution to the solar neutrino 
problem: the lack of success to date may merely reflect our lack of creativity. But if such 
a model exists - one that is consistent with our general observations of main-sequence 
evolution of solar-type stars and with helioseismology - it likely involves some new and 
subtle physics. 

3. THE DETECTION OF SOLAR NEUTRINOS 

Four solar neutrino experiments have now provided data, the Homestake ^^Cl ex- 
periment, the gallium experiments SAGE and GALLEX, and Kamiokande. The first 
three detectors are radiochemical, while Kamiokande records neutrino-electron elastic 
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scattering event- by-event. 

3.1 The Homestake Experiment 

Detection of neutrinos by the reaction ^^Cl(z/e,e)^^Ar was suggested independently by 
Pontecorvo (1946) and by Alvarez (1949). Davis's efforts to mount a 0.61 kiloton exper- 
iment using perchloroethylene (C2CI4) were greatly helped by Bahcall's demonstration 
(1964a, 1964b) that transitions to excited states in ^^Ar, particularly the superallowed 
transition to the analog state at 4.99 MeV, increased the *B cross section by a factor 
of 40. This suggested that Davis's detector would have the requisite sensitivity to de- 
tect neutrinos, thereby accurately determining the central temperature of the sun. 
The experiment (see Figure 7) was mounted in the Homestake Gold Mine, Lead, South 
Dakota, in a cavity constructed approximately 4850 feet underground [4900 meters wa- 
ter equivalent (m.w.e.)]. It has operated continuously since 1967 apart from a 17 month 
hiatus in 1985/86 caused by the failure of the circulation pumps. The result of 25 years 
of measurement is (Lande 1995) 

((70)37^ = 2.55 ±0.17 ±0.18 SNU (la) (6) 

which can be compared to the BP and TCL SSM predictions of 8.0 ± 1.0 SNU and 
6.4 ±1.4 SNU, respectively, all with la errors. As we will discuss below, the ^B and 
^Be contributions account for 77% (73%) and 15% (17%), respectively, of the BP (TCL) 
total. 

The experiment (Davis 1985; Davis 1993) depends on the special properties of ^^Ar: 
as a noble gas, it can be removed readily from perchloroethylene, while its half life (ti/2 = 
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35 days) allows both a reasonable exposure time and counting of the gas as it decays back 
to ^^Cl. Argon is removed from the tank by a helium purge, and the gas then circulated 
through a condensor, a molecular sieve, and a charcoal trap cooled to the temperature of 
hquid nitrogen. Typically ~ 95% of the argon in the tank is captured in the trap. (The 
efficiency is determined each run from the recovery results for a known amount of carrier 
gas, ^^Ar or ^^Ar, introduced into the tank at the start of the run.) When the extraction 
is completed, the trap is heated and swept by He. The extracted gas is passed through a 
hot titanium filter to remove reactive gases, and then other noble gases are separated by 
gas chromatography. The purified argon is loaded into a small proportional counter along 
with tritium-free methane, which serves as a counting gas. Since the electron capture 
decay of ^^Ar leads to the ground state of ^^Cl, the only signal for the decay is the 2.82 
keV Auger electron produced as the atomic electrons in ^^Cl adjust to fill the K-shell 
vacancy. The counting of the gas typically continues for about one year (~ 10 half lives). 

The measured cosmic ray-induced background in the Homestake detector is 0.06 ^^Ar 
atoms/ day while neutron-induced backgrounds are estimated to be below 0.03 atoms/day. 
A signal of 0.48 ± 0.04 atoms/day is attributed to solar neutrinos. When detector 
efficiencies, ^^Ar decays occurring in the tank, etc., are taken into account, the number 
of ^'^Ar atoms counted is about 25/year. 

A variety of careful tests of the argon recovery and counting efficiency have been 
made over the past 25 years. For example, known amounts of ^^Cl were introduced into 
the tank in order to check the recovery of ^^Ar, with a resulting yield of 100 ± 3%. It 
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has also been verified that ^^Ar produced in the tank by a fast neutron source, which 
induces (n,p) reactions followed by (p,n) on ^^Cl, is quantitatively recovered. However, 
the detector has never been calibrated directly with a neutrino source, despite studies of 
the feasibility of a ^^Zn source (Alvarez 1973). 

The significance of the Homestake results is due in part to an accurately determined 
^^Cl cross section. As the 814 keV threshold for exciting the ^^Ar ground state is above 
the pp endpoint, the detector is sensitive primarily to ^Be and neutrinos (see Table 
2). The cross section for '''Be neutrinos (and the weaker fluxes of pep and CNO cycle 
neutrinos) is determined by the known half life of ^^Ar. However, ^B neutrinos can 
generate transitions to many excited states below the particle breakup threshold in ^''Ar. 
The superallowed transition to the 4.99 MeV state, dominated by the Fermi matrix 
element of known strength, accounts for about 60% of the SSM cross section. The 
allowed transition strengths can be measured by observing the delayed protons following 
the P decay of ^^Ca, the isospin analog of the reaction ^'''Cl(fe,e) ^^Ar (Bahcall 1966). 
While it was believed that this measurement had been properly done many years ago 
(Poskanzer et al. 1966; Sextro, Gough, and Cerny 1974), the issue was not resolved 
until kinematically complete measurements were done recently (Garcia et al. 1991 and 
1995; Adelberger and Haxton 1987). The net result is a ^^Cl cross section believed to be 
accurate to about 3%. 
3.2 The Kamiokande Experiment 

The Kamiokande experiment (Hirata et al. 1988; Hirata et al. 1991) is a 4.5 kiloton 
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cylindrical imaging water Cerenkov detector originally designed for proton decay searches, 
but later reinstrumented to detect low energy neutrinos. It detects neutrinos by the 
Cerenkov light produced by recoihng electrons in the reaction 

i/^ + e ^ I/; + e'. (7) 

Both and heavy flavor neutrinos contribute, with cr(z/e)/cr(i/^) ~ 7. The inner volume 
of 2.14 kilotons is viewed by 948 Hamamatsu 20" photomultipher tubes (PMTs) pro- 
viding 20% photocathode coverage, and the surrounding 1.5m of water, serving as an 
anticounter, is viewed by 123 PMTs. The fiducial volume for solar neutrino measure- 
ments is the central 0.68 kilotons of water, the detector region most isolated from the 
high energy gamma rays generated in the surrounding rock walls of the Kamioka mine. 

In the conversion of the original proton decay detector to Kamiokande II, great effort 
was invested in reducing low energy backgrounds associated with radon and uranium. 
This included sealing the detector against radon inleakage and recirculating the water 
through ion exchange columns. The relatively shallow depth of the Kamioka mine (2700 
m.w.e.) leads to an appreciable flux of cosmic ray muons which, on interacting with ^^O, 
produce various short-lived spallation products. These [3 decay activities are vetoed by 
their correlation in time with the muons. The experimenters succeeded in lowering the 
detector threshold to 9 MeV and later to 7.5 MeV. Kamiokande III included improve- 
ments in the electronics and the installation of wavelength shifters around the PMTs to 
increase light collection and currently operates with a threshold of 7.0 MeV. 

Kamiokande II/III detects the high energy portion of the neutrino spectrum. Be- 
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tween December, 1985, and July, 1993, 1667 live detector days of data were accumulated. 
Under the assumption that the incident neutrinos are feS with an undistorted (3 decay 
spectrum, the combined Kamiokande II/III data set gives (Nakamura 1994) 

(p^^i^B) = (2.91 ± 0.08 ± 0.12) • (la) (8) 

corresponding to 51% of the BP and 63% of the TCL SSM predictions. The total number 
of detected solar neutrino events is 476^34. 

This experiment is remarkable in several respects. It is the first detector to measure 
solar neutrinos in real time. Essential to the experiment is the sharp peaking of the elec- 
tron angular distribution in the direction of the incident neutrino: this forward peaking, 
illustrated in Figure 8, allows the experimenters to separate solar neutrino events from 
an isotropic background. The unambiguous observation of a peak in the cross section 
correlated with the position of the sun is the first direct demonstration that the sun pro- 
duces neutrinos as a byproduct of fusion. Finally, although reaction (9) is a soft process 
where the recoil electron and scattered neutrino share the initial energy, the recoil elec- 
tron energy distribution provides some information on the incident neutrino spectrum. 
The recoil spectrum measured by Kamiokande II/IIl, shown in Figure 9, is consistent 
with an allowed incident neutrino spectrum, with the overall flux reduced as in Eq. 
(11). 

3.3 The SAGE and GALLEX Experiments 

Two radiochemical gallium experiments exploiting the reaction ^^Ga(i/e,e)^^Ge, SAGE 
and GALLEX, began solar neutrino measurements in January, 1990, and May, 1991, 
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respectively. SAGE operates in the Baksan Neutrino Observatory, under 4700 m.w.e. 
of shielding from Mount Andyrchi in the Caucasus, while GALLEX is housed in the 
Gran Sasso Laboratory at a depth of 3300 m.w.e. These experiments are sensitive 
primarily to the low-energy pp neutrinos, the flux of which is sharply constrained by 
the solar luminosity in any steady-state model of the sun (see Table 2). The gallium 
experiment was first suggested by Kuzmin (1966). In 1974 Ray Davis and collabora- 
tors began work to develop a practical experimental scheme. Their efforts, in which 
both GaCl3 solutions and Ga metal targets were explored, culminated with the 1.3- 
ton Brookhaven/Heidelberg/Rehovot/Princeton pilot experiment in 1980-82 that demon- 
strated the procedures later used by GALLEX (Hampel 1985). 

The primary obstacles to mounting the gallium experiments were the cost of the target 
and the greater complexity of the ^^Ge chemical extraction. The GALLEX experiment 
(Anselmann et al. 1992 and 1994a) employs 30 tons of Ga as a solution of GaCls in 
hydrochloric acid. After an exposure of about three weeks, the Ge is recovered as GeCU 
by bubbling nitrogen through the solution and then scrubbing the gas through a water 
absorber. The Ge is further concentrated and purified, and finally converted into GeH4 
which, when mixed with Xe, makes a good counting gas. The overall extraction efficiency 
is typically 99%. The GeH4 is inserted into miniaturized gas proportional counters, 
carefully designed for their radiopurity, and the Ge counted as it decays back to Ga 
(ti/2 = 11.43 d). As in the case of ^''Ar, the only signal for the Ge decay is the energy 
deposited by Auger electrons and x-rays that accompany the atomic rearrangement in 
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Ga. An important achievement of GALLEX has been the detection of both the K peak 
(10.4 keV) and L peak (1.2 keV). While 88% of the electron captures occur from the K 
shell, many of the subsequent K — > L x-rays escape the detector and some of the Auger 
electrons hit the detector walls. This produces a shift of the detected energy of events 
from the K to the L and M peaks. Thus the GALLEX L-peak counting capability almost 
doubles the ^^Ge detection efficiency. 

Gallium, like mercury, is a liquid metal at room temperature. SAGE (Abdurashitov 
1994; Bowles and Gavrin 1993) uses metallic gallium as a target, separating the ^^Ge by 
vigorously mixing into the gallium a mixture of hydrogen peroxide and dilute hydrochloric 
acid. This produces an emulsion, with the Ge migrating to the surface of the emulsion 
droplets where it is oxidized and dissolved by hydrochloric acid. The Ge is extracted 
as GeCl4, purified and concentrated, synthesized into GeH4, and further purified by 
gas chromatography. The overall efficiency, determined by introducing a Ge carrier, is 
typically 80%. The Ge counting proceeds as in GALLEX. SAGE began operations with 
30 tons of gallium, and now operates with 55 tons. The combined result for stage I (prior 
to September, 1992) and the first nine runs of stage II (9/92 - 6/93) is (Nico 1995) 

(a0)7iGe = 69+}} (stat)+5 (sys) SNU (la). (9) 

This result includes only events counted in the K peak. The counter and electronics 
improvements made at the start of stage II should permit L-peak events to be included, 
but no results have been announced as of December, 1994. The corresponding results for 
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GALLEX I (first 15 runs) and II (8/92-10/93) is (Anselmann et al. 1994) 



(a(/))nGe = 79 ± 10 ± 6 SNU (la). 



(10) 



GALLEX II solar neutrino runs continued until June, 1994, but results for the last portion 
of this period have not been reported. 

GALLEX II solar neutrino runs were interrupted in June, 1994, to permit an overall 
test of the detector with a ^^Cr source, which produces line sources of 746 keV (90%) 
and 426 keV (10%) neutrinos. The 1.67 MCi source was produced by irradiating ~ 36 
kg of chromium, enriched in ^'^Cr, in the Siloe reactor in Grenoble. Following exposure 
of the detector and recovery and counting of the produced ^^Ge, the ratio of measured 
^^Ge to expected was calculated (Anselmann et al. 1995), 



This is the first test of a solar neutrino detector with an terrestrial, low energy neutrino 
source. A similar source (~ 0.5 MCi) has been produced by the SAGE collaboration 
and was installed in their detector in December, 1994 (J.F. Wilkerson, private commu- 
nication). The higher Ga density of the SAGE detector increases the effectiveness of the 
source by about a factor of 2.5, helping to compensate for the weaker neutrino flux. 

The nuclear physics of the reaction ''^Ga(i/e,e)^^Ge is illustrated in Figure 10. As 
the threshold is 233 keV, the ground state and first excited state can be excited by pp 
neutrinos. However, as only those pp neutrinos within 12 keV of the endpoint can reach 
the excited state, the phase space for reaching this state is smaller by a factor of ~ 



= 1.04 ±0.12 (Icr). 



(11) 
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100. Thus the cross section is determined precisely by the measured electron capture 
lifetime of ^^Ge. In the BP and TCL calculations these neutrinos account for 54% and 
57% of the capture rate, respectively, each predicting 71 SNU. Because of this strong 
pp neutrino contribution, there exists a minimal astronomical counting rate of 79 SNU 
(Bahcall 1989) for the Ga detector that assumes only a steady-state sun and standard 
model weak interaction physics. This minimum value corresponds to a sun that produces 
the observed luminosity entirely through the ppl cycle. The rates found by SAGE and 
GALLEX are quite close to this bound. 

The ^Be neutrinos can excite the ground state and two excited states at 175 keV 
(5/2~) and 500 keV (3/2~). The SSM rates quoted above assume that these excited 
state transitions are much weaker than the ground state transition, contributing only 5% 
to the ^Be rate. The primary justification for this is the forward-angle (p,n) calibration 
of Gamow- Teller transitions to these states by Krofcheck et al. (1985). However, while 
the track record of (p,n) mappings of the broad features of the Gamow-Teller resonance 
is quite good, the technique is not generally considered a reliable test of the GT strength 
carried by an isolated state unless the transition to that state is quite strong (Austin et 
al. 1994). On the other hand, all three ^Be transitions are tested in the ^^Cr calibration, 
weighted by phase space factors quite similar to those of ^Be neutrinos (Haxton 1988a). 
Given other checks the GALLEX collaboration has made of the detector's overall chemical 
efficiency, another interpretation (Hata and Haxton 1995) of the source experiment is that 
it verifies that the excited states play a minor role in the ^Be capture rate. As the SSM 
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^Be capture rates of BP and TCL are both above 30 SNU, the SAGE/GALLEX results 
alone suggest some reduction in the low-energy pp and ^Be fluxes. 

The ^B neutrino capture rates (14 and 11 SNU in the BP and TCL SSMs, respectively) 
have been calculated from the GT profile deduced by Krofcheck et al. (1985). The 
corresponding total SSM rates for this detector are 132 SNU and 123 SNU (see Table 2). 
3.4 The Molvbdenum Experiment 

Recently an effort was made to do a solar neutrino experiment of a different kind, a 
geo chemical integration of the ^B neutrino flux over the past several million years (Cowan 
and Haxton 1982; Wolfsberg 1985). The reaction 

i/e Mo ^ e- Tc (ri/2 = 4.2 • lO^y) (12) 

occurring over geologic times can produce concentrations of '^'^Tc of ~ 10 atoms per 
gram of ^^Mo (abundance 24%). Because no stable isotopes of technetium exist, such 
a concentration is, in principle, measureable. Careful calculations of backgrounds from 
natural radioactivity and cosmic rays indicated they would be tolerable in a deeply buried 
ore body. 

The motivation of the experiment differs from others we have discussed: a compari- 
son of ^B flux averaged over an appreciable fraction of the Kelvin-Helmholtz time with 
contemporary measurements would test the long-term thermal stability of the core. 

The challenge is to isolate and quantitatively count ~ 10^ atoms of ^^Tc from 2600 
tons of raw ore (containing 13 tons of molybdenum) obtained from the one operating 
deep mine in North America, the Henderson Mine in Clear Creek County, Colorado. 
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Two major tasks faced the experimenters, isolating the Tc from this enormous quantity 
of ore and doing quantitative Tc mass spectrometry at the required level. The first goal 
appeared reachable as it is known that chemically analogous rhenium can be quantita- 
tively recovered from the gas stream of commercial molybdenum roasters. The second 
goal was reached after several years of work. 

The experiment was mounted by Los Alamos National Laboratory with the help 
of the AMAX Mining Corporation. The experimenters achieved a sensitivity to Tc at 
about the 100 SNU level, five times the expected SSM production rate. At that level 
a background appeared due to the commercial roaster's "memory" of recently roasted 
molybdenum from shallow mines in which cosmic-ray-induced Tc levels are high. While 
the plant memory decays with time, the necessary fiushing of the roaster with several 
weeks of Henderson ore was not practical commercially. The effort was abandoned in 
1988. 

It appears this experiment must await a factor of 10-100 further improvement in Tc 
mass spectrometry. At this sensitivity, the roasting of the M0S2 concentrate could be 
done under controlled laboratory conditions on a "table-top" scale of 10-100 kgs. 
4. PARTICLE PHYSICS SOLUTIONS 

In Section 2 it was shown that solar models which reduce the high energy neutrino 
fiux tend to enhance the ^Be/^B fiux ratio, contradicting the results of the Homestake, 
SAGE/GALLEX, and Kamiokande experiments. 

If the source of the solar neutrino problem is not solar, the remaining possibilities are 
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experimental error or nonstandard particle physics. Several researchers have recently con- 
sidered the consequences of ignoring one of the three experiments just discussed (Parke 
1995; Bahcall 1994; Kwong and Rosen, 1994; Bahcall and Bethe 1990). Substantial 
discrepancies between SSM predictions and the remaining experiments persist. For ex- 
ample, Figure 11 shows Parke's results when only the Kamiokande and SAGE/GALLEX 
constraints are retained: a discrepancy of almost 3-4 a remains, depending on the choice 
of SSM. The corresponding results for other combinations show even larger inconsisten- 
cies. If two experiments must be flawed to account for the discrepancy with the SSM, 
this scenario becomes somewhat less credible. 

The second alternative, physics beyond the standard model of electroweak inter- 
actions, would have the most far-reaching consequences. Particle physics solutions of 
the solar neutrino problem include neutrino oscillations, neutrino decay, neutrino mag- 
netic moments, and weakly interacting massive particles. Among these, the Mikheyev- 
Smirnov-Wolfenstein effect — neutrino oscillations enhanced by matter interactions — is 
widely regarded as the most plausible. 
4.1 Neutrino Masses and Vacuum Oscillations 

One odd feature of particle physics is that neutrinos, which are not required by any 
symmetry to be massless, nevertheless must be much lighter than any of the other known 
fermions. For instance, the current limit on the mass is ^ 5 eV. The standard model 
requires neutrinos to be massless, but the reasons are not fundamental. Dirac mass terms 
mo, analogous to the mass terms for other fermions, cannot be constructed because the 
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model contains no right-handed neutrino fields. Neutrinos can also have Majorana mass 
terms 

^i'rnLJ^L and u^ruRi^R. (13) 

where the subscripts L and R denote left- and right-handed projections of the neu- 
trino field u, and the superscript c denotes charge conjugation. The first term above is 
constructed from left-handed fields, but can only arise as a nonrenormalizable effective 
interaction when one is constrained to generate rriL with the doublet scalar field of the 
standard model. The second term is absent from the standard model because there are 
no right-handed neutrino fields. 

None of these standard model arguments carries over to the more general, unified the- 
ories that theorists believe will supplant the standard model. In the enlarged multiplets 
of extended models it is natural to characterize the fermions of a single family, e.g., fe, 
e, u, d, by the same mass scale mij. Small neutrino masses are then frequently explained 
as a result of the Majorana neutrino masses. In the seesaw mechanism. 



tud 



y niR J 



(14) 



2 

Diagonalization of the mass matrix produces one light neutrino, mught ~ and one 



unobservably heavy, mheavy ~ ^r- The factor {mD/mji) is the needed small parameter 
that accounts for the distinct scale of neutrino masses. The masses for the 1/^, v^, and 
Ut- are then related to the squares of the corresponding quark masses m^, rric, and m^. 
Taking niR ~ 10^^ GeV, a typical grand unification scale for models built on groups like 
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so (10), the seesaw mechanism gives the crude relation 

m^^ : m^^ : m^^ ^ 2 ■ 10"^^ : 2 • 10"^ : 3 • lO'^eV. (15) 

The fact that solar neutrino experiments can probe small neutrino masses, and thus pro- 
vide insight into possible new mass scales thr that are far beyond the reach of direct 
accelerator measurements, has been an important theme of the field (Babu and Mo- 
hapatra 1993; Bludman, Kennedy, and Langacker, 1992; Dimopoulos, Hall, and Raby 
1993). 

Another expectation is that neutrinos from the different families mix, just as quark 
mixing is observed in hyperon and nucleon /3 decays. If we consider the two-flavor case 
for simplicity, the mass eigenstates and I1/2) (with masses mi and 777.2) are related to 
the weak interaction eigenstates by 

l^e) — cos^.y|i/i) -|- sin^„|z/2) 

= — sin6'^|i/i) -h cos6'^|i/2) (16) 

where 9y is the (vacuum) mixing angle. The two mass eigenstates comprising the then 
propagate with different phases in vacuum, leading to flavor oscillations. The probability 
that a will remain a Ug after propagating a distance x is 

P^^{x)^l-sin'29,sm'(^-^] ^ (17) 

where E is the neutrino energy and 6m^ = m\ — ml. (When one properly describes the 
neutrino state as a wave packet, the large-distance behavior follows from the eventual 
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separation of the mass eigenstates.) If the the oscillation length 

is comparable to or shorter than one astronomical unit, a reduction in the flux would 
be expected in terrestrial neutrino oscillations. The suggestion that the solar neutrino 
problem could be explained by neutrino oscillations was first made by Pontecorvo (1958), 
who pointed out the analogy with Kq Kq oscillations. From the point of view of 
particle physics, the sun is a marvelous neutrino source. The neutrinos travel a long 
distance and have low energies (~ 1 MeV), implying a sensitivity to 

In the seesaw mechanism, ~ m|, so neutrino masses as low as ~ l^'^eV could 
be probed. In contrast, terrestrial oscillation experiments with accelerator or reactor 
neutrinos are typically limited to brv? > CleV^^. 

Prom Eq. (17) one expects vacuum oscillations to affect all neutrino species equally, 
if the oscillation length is small compared to an astronomical unit. This appears to 
contradict observation, as the pp flux may not be significantly reduced. Furthermore, 
the theoretical prejudice that should be small makes this an unlikely explanation of 
the significant discrepancies with SSM ^Be and fiux predictions. 

The first objection, however, can be circumvented in the case of "just so" oscillations 
where the oscillation length is comparable to one astronomical unit (Glashow and Krauss 
1987). In this case the oscillation probability becomes sharply energy dependent, and one 
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can choose SmP' to preferentially suppress one component (e.g., the monochromatic ^Be 
neutrinos) . This scenario has been explored by several groups and remains an interesting 
possibility. However, the requirement of large mixing angles remains. 
4.2 The Mikhevev-Smirnov-Wolfenstein Mechanism 

The community's view of neutrino oscillations changed radically when Mikheyev and 
Smirnov (1985 and 1986) showed that the density dependence of the neutrino effective 
mass, a phenomenon first discussed by Wolfenstein (1978), could greatly enhance oscil- 
lation probabilities: a is adiabatically transformed into a as it traverses a critical 
density within the sun. It became clear that the sun was not only an excellent neutrino 
source, but also a natural regenerator for cleverly enhancing the effects of flavor mixing. 

While the original work of Mikheyev and Smirnov was numerical, their phenomenon 
was soon understood analytically as a level-crossing problem. If one writes the neutrino 
wave function in matter as 

\u{x)) = ae{x)\u^) + a^(x)|i/^) 

where x is the coordinate along the neutrino's path, the evolution of ae{x) and a^(x) is 
governed by 

/ de 
\ai_i 

where G^^ is the weak coupling constant and p{x) the solar electron density. If p{x) — 0, 
Eq. (20) can be trivially integrated to give the vacuum oscillation solution (Eq. (17)). 
The new contribution to the diagonal elements, 2E\/2Gfp{x), represents the effective 
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'p{x) — Sm^ cos 26v Sm^ sin 29v \ 

h - 2E^/2Gfp{,x) + cos 2^^; / 



(20) 



contribution to that arises from neutrino-electron scattering. The indices of refraction 
of electron and muon neutrinos differ because the former scatter by charged and neutral 
currents, while the latter have only neutral current interactions. The difference in the 
forward scattering amplitudes determines the density-dependent splitting of the diagonal 
elements of Eq. (20). 

It is helpful to rewrite Eq.(20) in a basis consisting of the hght and heavy local mass 
eigenstates (i.e., the states that diagonalize the right-hand side of Eq. (20)), 

luLix)) = cos9{x)\i/e) - sin6l(x)|i/^) 



Wh{x)) = sin6'(a;)|z/e) -|- cos6'(a;)|z/^). 



(21) 



The local mixing angle is defined by 

sin 29 (x) 

cos2^(x) 



sin 29v 



^Jx\x)+shi^29v 

-Xix) 



(22) 



X\x) + sin^ 29v 

where X{x) = 2y/2GFp{x)E /Sm^ - cos29v. Thus 9{x) ranges from 9v to 7r/2 as the 



density p{x) goes from to oo. 
If we define 



|i/(x)) = aH{x)\vH{x)) aL{x)\vL{x)), 



Eq. (20) becomes 
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(23) 
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with the sphtting of the local mass eigenstates determined by 



2A(a;) 



~2E 



(24) 



and with mixing of these eigenstates governed by the density gradient 



Note that the splitting achieves its minimum value, ^|rsin2^^, at a critical density 



(25) 



Pc = p{Xc) 



cos 29., 



(26) 



that defines the point where the diagonal elements of Eq. (20) cross. 

Equation (23) can be trivially integrated if the splitting of the diagonal elements is 
large compared to the off-diagonal elements, 



7(x) 



a[Xj 



sin' 2ev 5m^ 1 [X{xf + sin' 29^] 

' T~r~\ : o — 1, 



cos2^^; 2E \±Me1\ sm^29„ 

I Pc dx I 



a condition that becomes particularly stringent near the crossing point. 



7c = l{Xc) 



sin'26't, Sm^ 



cos 26*,, 2E 



1 dp{x) I 
pc dx I ^=^c 



> 1. 



(27) 



(28) 



The resulting adiabatic electron neutrino survival probability, valid when 7c ^ 1, is 



padiab ^ 1 ^ 1 (,os 20 cos 29i 

I/e 2 2 



(29) 



where 9i — 9{xi) is the local mixing angle at the density where the neutrino was produced. 
Eq. (29) was first discussed by Bethe (1986) (also see Messiah 1986). 



36 



The physical picture behind this derivation is illustrated in Figure 12. One makes the 
usual assumption that, in vacuum, the is almost identical to the light mass eigenstate, 
1/1,(0), i.e., nil < and cos^^ ~ 1. But as the density increases, the matter effects 
make the heavier than the i/^, with — > VHi^) as p{x) becomes large. The special 
property of the sun is that it produces v^s at high density that then propagate to the 
vacuum where they are measured. The adiabatic approximation tells us that if initially 

~ ^h{x), the neutrino will remain on the heavy mass trajectory provided the density 
changes slowly. That is, if the solar density gradient is sufficiently gentle, the neutrino 
will emerge from the sun as the heavy vacuum eigenstate, yui^) ~ ^n- This guarantees 
nearly complete conversion of i/gS into i/^s, producing a flux that cannot be detected by 
the Homestake or SAGE/GALLEX detectors. 

But this does not explain the curious pattern of partial flux suppressions of Eq. 
(2). The key to this is the behavior when 7c < 1. Eqs. (28) and (29) show that 
the critical region for nonadiabatic behavior occurs in a narrow region (for small 9^) 
surrounding the crossing point, and that this behavior is controlled by the derivative of 
the density. This suggests an analytic strategy for handling nonadiabatic crossings: one 
can replace the true solar density by a simpler (integrable!) two-parameter form that 
is constrained to reproduce the true density and its derivative at the crossing point Xc- 
Two convenient choices are the linear {p{x) — a -\- hx) and exponential {p{x) — ae"^^) 
profiles. As the density derivative at Xc governs the nonadiabatic behavior, this procedure 
should provide an accurate description of the hopping probability between the local mass 
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eigenstates when the neutrino traverses the crossing point. The initial and ending points 
Xi and Xf for the artificial profile are then chosen so that p{xi) is the density where the 
neutrino was produced in the solar core and p{xf) = (the solar surface), as illustrated 
in Fig. 13. Since the adiabatic result (Eq. (29)) depends only on the local mixing 
angles at these points, this choice builds in that limit. Eq. (20) can then be integrated 
exactly for linear and exponential profiles, with the results given in terms of parabolic 
cylinder and Whittaker functions, respectively. This treatment, called the finite Landau- 
Zener approximation (Haxton 1987; Petcov 1988) has been used extensively in numerical 
calculations. 

We derive a simpler ("infinite") Landau-Zener approximation (Landau 1932; Zener 
1932) by observing that the nonadiabatic region is generally confined to a narrow region 
around Xc, away from the endpoints Xi and Xf. We can then extend the artificial profile 
to X — ±oo, as illustrated by the dashed lines in Fig. 13. As the neutrino propagates 
adiabatically in the unphysical region x < Xi, the exact soluation in the physical region 
can be recovered by choosing the initial boundary conditions 

aii—oo) — — a^(— oo) = cos^je~*-^-°°'^*-^'*'^^ 

ani—oo) — ae(— oo) = sin^jc'-^-'^'^^^^'^^ (30) 

That is |i/(— oo)) will then adiabatically evolve to li'ixi)) — li^e) as x goes from — oo to 
Xi- The unphysical region x > Xf can be handled similarly. 

With some algebra a simple generalization of the adiabatic result emerges that is valid 
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for all dm^/E and 9^ 



^ + ^ cos 2^, cos2^i(l - 2Phop) (31) 



where Phop is the probability of hopping from the heavy mass trajectory to the light 
trajectory on traversing the crossing point. For the linear approximation to the density, 

P^^ = e-^^^^/^ (32) 

As it must by our construction, P,^^ reduces to Y"^^^^ for 7c ^ 1. The linear Landau- 
Zener asymptotic hopping probability Ph"p = e~'^'^''l'^ was derived by Haxton (1986) and 
independently by Parke (1986), who married this approximation to the adiabatic one to 
get Eq. (31). The exponential probability was first obtained by Petcov (1988), 



g— 7r<5(l— cos 20^) g- 



■27r<5 



where 5 = ^^^§7- Note that for small Q^, 5(1 - cos2e^) ^ 7e/2 and 5 ^ 7c/(2^^)^ so 
that Eqs. (32) and (33) then coincide. When the crossing becomes nonadiabatic (e.g., 
7c -C 1 in Eq. (32)), the hopping probability goes to 1, allowing the neutrino to exit the 
sun on the light mass trajectory as a i/g, i.e., no conversion occurs. 

Thus there are two conditions for strong conversion of solar neutrinos: there must be 
a level crossing (that is, the solar core density must be sufficient to render ~ T^Hi^i) 
when it is first produced) and the crossing must be adiabatic. The first condition requires 
that Sm^/E not be too large, and the second 7c ^ 1- The combination of these two 
constraints, illustrated in Figure 14, defines a triangle of interesting parameters in the 

39 



^ — sin^ 29y plane, as Mikheyev and Smirnov and others (Rosen and Gelb 1986) found 
by numerically integrating Eq. (20). A remarkable feature of this triangle is that strong 
i^e t^fi conversion can occur for very small mixing angles (sin^ 2^ ~ 10~^), unlike the 
vacuum case. 

One can envision superimposing on Fig. 14 the spectrum of solar neutrinos, plotted 
as a function of for some choice of Sm"^. Since Davis sees some solar neutrinos, the 
solutions must correspond to the boundaries of the triangle in Figure 14. The horizontal 
boundary indicates the maximum for which the sun's central density is sufficient 
to cause a level crossing. If a spectrum properly straddles this boundary, we obtain a 
result consistent with the Homestake experiment in which low energy neutrinos (large 
1/E) lie above the level-crossing boundary (and thus remain t'e's), but the high-energy 
neutrinos (small 1/E) fall within the unshaded region where strong conversion takes 
place. Thus such a solution would mimic nonstandard solar models in that only the 

neutrino flux would be strongly suppressed. The diagonal boundary separates the 
adiabatic and nonadiabatic regions. If the spectrum straddles this boundary, we obtain 
a second solution in which low energy neutrinos lie within the conversion region, but the 
high-energy neutrinos (small 1/E) lie below the conversion region and are characterized 
by 7 -C 1 at the crossing density. (Of course, the boundary is not a sharp one, but 
is characterized by the exponential of Eq. (32)). Such a nonadiabatic solution is quite 
distinctive since the flux of pp neutrinos, which is strongly constrained in the standard 
solar model and in any steady-state nonstandard model by the solar luminosity, would 
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now be sharply reduced. Finally, one can imagine "hybrid" solutions where the spectrum 
straddles both the level-crossing (horizontal) boundary and the adiabaticity (diagonal) 
boundary for small d, thereby reducing the ^Be neutrino flux more than either the pp or 
fluxes. 

What are the results of a careful search for MSW solutions satisfying the Homestake, 
Kamiokande, and SAGE/GALLEX constraints? Figure 15 is a calculation by Hata (1995) 
(also see Hata and Langacker 1994) for flavor oscillations that includes the effects of 
terrestrial regeneration. (MSW effects can occur as the neutrinos pass through the earth.) 
The preferred (in the sense of minimizing the x^) solution, corresponding to a region 
surrounding drri^ ~ 6 • lO'^eV"^ and sin^ 29y ~ 6 • 10~^, is the hybrid case described 
above. It is commonly called the small-angle solution. A second, large-angle solution 
exists, corresponding to Sm^ ~ 10~^eV^ and sin^ 29^ ~ 0.6, but this region of Figure 15 
has shrunk as the precision of the gallium experiments improve. 

These solutions can be distinguished by their characteristic distortions of the solar 
neutrino spectrum. The survival probabilities P^^^{E) for the small- and large-angle 
parameters given above are shown as a function of E in Figure 16. 

The calculations of Figure 14 assume flavor oscillations into a or i/^. This influences 
the interpretation of the Kamiokande experiment, as heavy flavor neutrinos contribute to 
elastic scattering. Another possibihty is an oscillation into a sterile neutrino. The large- 
angle solution is then ruled out by the Kamiokande requirement of a large Ue survival 
probability (see, for example, Hata 1995; Barger, Deshpande, Pal, PhiUips, and Whis- 
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nant, 1991). It is also ruled by the neutrino counting limit from big bang nucleosynthesis 
(Barbieri and Dolgov 1991; Enqvist el al. 1990; Shi, Schramm, and Fields 1993). 

The MSW mechanism provides a natural explanation for the pattern of observed 
solar neutrino fluxes. While it requires profound new physics, both massive neutrinos 
and neutrino mixing are expected in extended models. The preferred solutions correspond 
to ~ 10~^ eV^, and thus are consistent with 777.2 ~ few -lO"^ eV. This is a typical 
Vj- mass in models where 777i^ ~ ?77gut- On the other hand, if it is the participating 
in the oscillation, this gives ttir ~ 10^^ GeV and predicts a heavy v^. ~ 10 eV (Bludman, 
Kennedy, and Langacker 1992). Such a mass is of great interest cosmologically as it 
would have consequences for supernova neutrinos (Fuller et al. 1992; Qian et al. 1993), 
the dark matter problem, and the formation of large-scale structure. 

If the MSW mechanism proves not to be the solution of the solar neutrino problem, 
it still will have greatly enhanced the importance of solar neutrino physics: the existing 
experiments have ruled out large regions in the Sm^ — sin^ 29^ plane (corresponding 
to nearly complete Uf. — > i/^ conversion) that remain hopelessly beyond the reach of 
accelerator neutrino oscillation experiments. 
4.3 Other Particle Phvsics Scenarios 

Several other intriguing particle physics phenomena could affect the solar neutrino 
puzzle. The upper bound established in the earliest ^^Cl runs was consistent with the 
complete absence of solar neutrinos, prompting the suggestion that the i/g might decay 
before reaching earth (Bahcall, Cabibbo, and Yahil 1972). This requires a neutrino mass 
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and a sufficiently fast decay mode. The suggestion is less appealing given the present 
pattern of fluxes: if the life time is arranged to allow some neutrinos to survive, one 
expects the low-energy neutrino flux to be more severely suppressed than the high-energy 
neutrinos. There is also the constraint from supernova SN1987A, where PgS successfully 
traveled 50 kpc. However, one can still wiggle out of both objections if the neutrino 
decay is catalyzed by matter effects 

^^Msw ~^ decay products (34) 

as Raghavan, He, and Pakvasa suggested (1988). The spectrum distortions do not neces- 
sarily mimic the MSW mechanism since the decay probability can depend on the neutrino 
energy. 

Weakly interacting massive particles (WIMPs) were suggested (Faulkner and Gilliland 
1985; Spergel and Press 1985) as a simultaneous solution to the solar neutrino and dark 
matter problems. If a heavy neutral particle has a mean free path (at solar densities) 
comparable to the solar diameter, it can lose energy in transit and be captured in the sun's 
gravitational field. Once the sun accumulates a sufficiently dense cloud of such particles, 
they contribute to energy transport in the sun, thus violating the SSM assumption that 
only radiative transport is important in the core. The WIMP can pick up energy by 
scattering off faster core nucleons, then lose it by rescattering in the cooler outer layers 
of the sun. The additional transport lowers the core temperature. While this suggestion is 
clever, a simple lowering of is no longer sufficient to reconcile the SSM with experiment. 

Perhaps the most interesting possibility, apart from the MSW mechanism, was stimu- 
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lated by suggestions that the ^^Cl signal might be varying with a period comparable to the 
11-year solar cycle. While the evidence for this has weakened, the original claims gener- 
ated renewed interest in neutrino magnetic moment interactions with the solar magnetic 
field. 

The original suggestions by Cisneros (1971) and by Okun, Voloshyn, and Vysotsky 
(1986) envisioned the rotation 

^ea (35) 

producing a right-handed neutrino with sterile interactions in the standard model. With 
the discovery of the MSW mechanism, it was realized that matter effects would break 
the vacuum degeneracy of the ^cr, suppressing the spin precession. Lim and 

Marciano (1988) pointed out that this difficulty was naturally circumvented for 

(36) 

as the different matter interactions of the and can compensate for the vacuum 
i^e — Vfj, mass difference, producing a crossing similar to the usual MSW mechanism. Such 
spin-flavor precession can then occur at full strength due to an off-diagonal (in flavor) 
magnetic moment. 

There has been a great deal of clever work on this problem (Minakata and Nunokawa 
1989; Balantekin, Hatchell, and Loreti 1990). A very strong limit on both diagonal and 
off-diagonal magnetic moments is imposed by studies of the red giant cooling process of 
plasmon decay into neutrinos 

7* ^ PiPj. (37) 
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The result is \fj,ij\ <^ 3 • 10~^^/xb, where fj,B is an electron Bohr magneton (RafFaelt 1990). 
With this bound, solar magnetic field strengths of Bq ^ lO^G are needed to produce 
interesting effects. Since the location of the Lim-Marciano level crossing depends on the 
neutrino energy, such fields have to be extensive to affect an appreciable fraction of the 
neutrino spectrum. It is unclear whether these conditions can occur in the sun. 

There are interesting, related phenomena involving the effects of solar density fiuc- 
tuations or currents on the MSW mechanism: one can drive i/g — > i/^ oscillations in the 
absence of a level crossing by harmonic density perturbations (not unlike adiabatic fast 
passage in nuclear magnetic resonance experiments) (Haxton and Zhang 1991; Schafer 
and Koonin 1987; Krastev and Smirnov 1989). The effects of "white noise" density fluctu- 
ations on the MSW mechanism have also been examined recently (Loreti and Balantekin 
1994). Such fluctuations generate a flavor analog of stochastic spin depolarization, a 
phenomenon familiar in atomic physics. 
5. NEW EXPERIMENTS 

The MSW mechanism has had a particularly strong impact because it was discov- 
ered at a time when new data (SAGE/GALLEX, Kamiokande, helioseismology) were 
eliminating many competing solutions to the solar neutrino problem. The physics of the 
MSW mechanism is both simple and elegant, which accounts for much of its appeal. But 
the most important attribute of this solution is that it can be definitively tested. The 
favored small-angle solution produces a distinctive distortion in the solar neutrino spec- 
trum. Furthermore, if the oscillation is into another fiavor (rather than a sterile state). 
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the missing neutrinos can be found through their neutral current interactions. These tests 
will be made by two high-statistics, direct-counting directors now under construction. 
5.1 Tlic Sudbury Neutrino Observatory 

A water Cerenkov detector of a different type is under construction deep (5900 m.w.e.) 
within the Creighton ^^9 nickel mine at Sudbury, Ontario, Canada (Ewan et al. 1987; 
Aardsma et al. 1987; Chen 1985). The central portion of the detector is an acrylic vessel 
containing 1 kiloton of heavy water, D2O. This is surrounded by five meters of light water 
to protect the inner detector from neutrons and gammas. The detector is viewed by 9500 
20-cm PMTs, providing 56% photocathode coverage (Figure 17). 

The D2O introduces two new channels. The charged current breakup reaction 

u^ + D^p + p + e- (38) 

produces a recoil electron which carries off almost all of the final-state kinetic energy. As 
the Gamow- Teller strength is concentrated very close to the p-|-p threshold, 1.44 MeV, 
the electron and neutrino energies are related by + 1.44 MeV. Thus, neutrino 

spectrum distortions should show up clearly in the measured electron energy distribution. 
As the GT strength in the deuteron is equivalent to about one-third of a free neutron, 
the anticipated counting rates are high. For an electron detection threshold of 5 MeV 
and a neutrino flux equal to 51% of the BP SSM value, 3300 events will be recorded 
each year. 

A second channel is sensitive equally to neutrinos of any flavor, 

u^ + D^u'^ + n + p (39) 
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and thus will be crucial in testing whether flavor oscillations have occurred. The antici- 
pated event rate is approximately 2000/year in the BP SSM. The addition of MgCl2 to 
the D2O at a concentration of 0.2-0.3% allows the neutrons to be observed by ^^Cl(n,7). 
The Cerenkov light produced by the showering of the 8.6 MeV capture 7 ray will add to 
the signal from the charged current reaction. By operating the detector with and without 
salt, the experimenters will separate the charged and neutral current signals. The SNO 
collaboration also plans to deploy proportional counters filled with ^He to exploit the 
neutron-specific charge-exchange reaction ^He(n,p)'^H. With such detectors, SNO will be 
sensitive to neutral current events at all times. 

The detection of ~ 8 neutrons/day in a kiloton detector places extraordinary con- 
straints on radiopurity. For example, a potentially serious background source is the 
photodisintegration of deuterium by energetic photons from U and Th chains. The ex- 
perimental goal is concentrations of <^ 10~^^ grams of U and Th per gram of D2O. 

SNO is scheduled to begin operations in mid-1996. 
5.2 Super kamiokande 

Super kamiokande will be a greatly enlarged version of Kamiokande II/III with im- 
proved threshold (5 MeV) and energy and position resolution (Totsuka 1987 and 1990; 
Takita 1993). It is currently under construction in the Kamioka mine at a depth of 2700 
m.w.e. (See Figure 18). 

The fiducial volume for detecting solar neutrinos will be 22 kilotons, compared to 
0.68 kilotons in the existing detector. This plus the improved threshold will increase 
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the detection rate for neutrino-electron scattering by a factor of ~ 90, to 8400/year. 
Despite the soft kinematics of the Vf. — e reaction, the experimenters beheve the high 
statistics will allow them to distinguish the spectral distortions produced by competing 
MSW solutions. 

Because elastic scattering is sensitive to both Vf. and heavy-flavor neutrinos (in the ra- 
tio of 7:1), an accurate SNO determination of the v^. spectrum will allow Super kamiokande 
to extract the spectrum of or v^%. 

Superkamiokande construction is scheduled to be completed in March, 1996. 
5.3 Other Future Detectors 

The Borexino collaboration (Raghavan 1991; Campanella 1992) has proposed a 0.3 
kiloton liquid scintillator for installation in the Gran Sasso Laboratory. The experi- 
menters hope to detect ^Be neutrinos by — e scattering. The detection of very low 
energy recoil electrons places stringent constraints on U, Th, K, and other activities in 
the detector, e.g., ^ 10~^^ g U/g, including a requirement for continuous purification. 
The experimenters will evaluate background problems in a test facility now under con- 
struction, and scheduled to be completed by the end of 1995. The anticipated counting 
rate for the full-scale detector is ~ 18,000 ^Be neutrino events/year for the BP SSM. 

A high-counting-rate twin of the ^^Cl detector utilizing the reaction ^^^I(i/e,e)^^^Xe 
has been funded recently and is under construction in the Homestake mine (Lande 1993; 
Haxton 1988b). With a threshold of 664 keV, the detector is primarily sensitive to ''Be 
and ^B neutrinos. The initial Homestake detector will contain 100 tons of iodine as 
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a solution of Nal. A smaller version of this detector was recently used at the LAMPF 
beamstop to measure the ^^^I cross section for stopped muon decay VeS. Calibration of 
the ^Be neutrino cross section by an ^^Ar neutrino source is also planned. The 100-ton 
detector is scheduled to be completed by the end of 1995, with plans for an expansion to 
1 kiloton afterwards. 

A 5-kiloton liquid argon time projection chamber, ICARUS II, has been proposed for 
Gran Sasso (Rubia 1985). In addition to u — e scattering, the charged current channel 



will allow the experimenters to measure the shape of the high-energy portion of the 
spectrum. 

There are a number of important development efforts underway that focus on new 
technologies for the next generation of solar neutrino detectors. The reader is referred 
to the recent review by Lanou (1995). There are significant challenges motivating these 
efforts, e.g., neutrino detection by coherent scattering off nuclei and real-time detectors 
for pp neutrinos, such as the superfluid '^He detector HERON (Handler et al. 1992) and 
the high-pressure helium time projection chamber HELLAZ (Laurenti et al. 1994). 



The successes of the Kamiokande and SAGE/GALLEX experiments have led to a 
more complicated solar neutrino problem. The apparent strong suppression of the ^Be 

flux (negative in most unconstrained fits!) is not a result expected by those who favored 
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6. OUTLOOK 
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nonstandard solar model solutions to the ^"Cl puzzle. Perhaps this is telling us that solar 
modelers have not been sufficiently inventive in modifying the SSM. But it may also be 
a push to look elsewhere for the solution. 

We do have one candidate solution that works extremely well, the MSW mechanism. 
The required new physics has deep implications for particle physics and cosmology. Yet 
this physics is not exotic - the requirements of massive neutrinos and mixing are com- 
mon assumptions in extensions of the standard electroweak model. The elegance of this 
solution makes it difficult to maintain one's scientific skepticism: the notion that the sun 
was perfectly designed to enhance the mixing of neutrinos with GUT-scale seesaw masses 
has too great an emotional appeal. 

Fortunately the solution to the solar neutrino problem does not have to be decided 
by community vote; the issue will be resolved by hard-nosed experimentation. SNO and 
Superkamiokande are just a year away, and they may crack this 30-year-old problem. 
Yet these are difficult experiments, and the physics they address is fundamental to two 
of our standard models (particle and solar). It is prudent for the community to continue 
to seek cross checks on these and other future measurements. As has proven true in the 
past, there is no guarantee that the SNO and Superkamiokande results will conform to 
our expectations. 

I am indebted to the participants of the Seattle Solar Modehng Workshop for much 
of the background material that was incorporated into this paper. It is a special pleasure 
to thank J. Bahcall, A. B. Balantekin, B. Cleveland, R. Davis Jr., S. Degl'Innocenti, 
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N. Hata, C. Johnson, K. Lande, K. Nakamura, S. Parke, P. Parker, Y.-Z. Qian, H. 
Robertson, Y. Totsuka, and J. Wilkerson for helpful discussions and assistance. This 
work was supported in part by the U.S. Department of Energy and by NASA under 
grant #NAGW2523. 
Note added in proof (3/24/95) 

Six of the eight planned exposures have been completed in the SAGE ^^Cr source 
experiment. The extractions have proceeded routinely, and an announcement of the 
results is expected later this year. 

The BP SSM has been updated by the inclusion of both hehum and heavy element 
diffusion (Bahcall and Pinsonneault 1995). The resulting ®B neutrino flux is 6.62E6 
cm~^s~^, an increase of 16% (31%) relative to the BP results with He diffusion (without 
diffusion) given in Table 1. The "best value" ^^Cl and ^^Ga capture rates are 9.3 and 137 
SNU, respectively. The present-day helium surface value is 0.247, in good agreement with 
the hehoseismology value of 0.242. The depth of the convective zone is also in excellent 
agreement with the value deduced from p-mode oscillation data. 
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Figure Captions 

Figure 1: The solar pp chain. 

Figure 2: The flux densities (sohd hues) of the principal P decay sources of solar neutrinos 
of the standard solar model. The total fluxes are those of the BP SSM. The ''Be and pep 
electron capture neutrino fluxes (dashed lines) are discrete and given in units of cm~^s"~^. 

Figure 3: The ^Be(p, 7)^B S-factors as measured by Kavanagh et al. (1969) and by 
Filippone et al. (1983). For each data set, two theoretical extrapolations, reflecting 
different choices for the strong potentials, to S(0) are shown (Johnson, Kolbe, Koonin, 
and Langanke 1992). The enlargement of the error bars is a correction by these authors 
to account for the systematic differences in the two data sets. 

Figure 4: The dots represent the ^Be and ^B fluxes resulting from the 1000 SSMs of 
Bahcall and Ulrich (1988), with smaller SSM uncertainties added as in Bahcall and 
Haxton (1989). The 90 and 99% c.l. error eUipses are shown. 

Figure 5: The response of the pp. Be, and ^B fluxes to the indicated variations in solar 
model input parameters, displayed as a function of the resulting central temperature Tg. 
From Castellani, Degl'Innocenti, Fiorentini, Lissia, and Ricci (1994). 

Figure 6: The fluxes allowed by the combined results of the Homestake, SAGE/GALLEX, 
and Kamiokande experiments are compared to the results of SSM variations and various 
nonstandard models. The sohd line is the Tc power law of Eq. (3). From Hata (1995). 
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Figure 7: The Homestake ^ CI solar neutrino experiment. This schematic is from the 
Davis, Harmer, and Hoffman (1968) paper reporting the first results from their experi- 
ment. 

Figure 8: Angular distribution of recoil electrons from Kamiokande II and III showing the 
excess at forward angles that is attributed to solar neutrinos. Electrons with apparent 
energies between 7 and 20 MeV are included. The upper histogram is the SSM prediction 
of Bahcall and Ulrich (1988) superimposed on an isotropic background, while the lower 
histogram is the best fit. From K. Nakamura (1994). 

Figure 9: The Kamiokande II/III recoil electron energy spectrum compared to the 
SSM prediction (solid histogram) and to the SSM prediction rescaled by a factor of 
0.5 (dashed). The last bin corresponds to electron apparent energies between 14 and 20 
MeV. From K. Nakamura (1994). 

Figure 10: Level scheme for ^^Ge showing the excited states that contribute to absorption 
of pp, '^Be, ^^Cr, and neutrinos. The ™Ge + n break-up threshold is 7.4 MeV. 

Figure 11: The ^Be and ^B fluxes determined by the SAGE/GALLEX and Kamiokande 
experiments are compared to the predictions of the BP and TCL SSMs. From Parke 
(1995). 

Figure 12: Schematic illustration of the MSW level crossing. The dashed hues correspond 
to the electron-electron and muon-muon diagonal elements of the matrix in Eq. (20). 
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Their intersection defines the level crossing density pc- The sohd lines are the trajectories 
of the light and heavy local mass eigenstates. If the electron is produced deep in the solar 
core and propagates adiabatically, it will follow the heavy mass trajectory, emerging from 
the sun as a i/^. 

Figure 13: The top figure illustrates, for one choice of sin^^„ and 5m^, that the region 
of nonadiabatic propagation (solid line) is usually confined to a narrow region about the 
crossing point Xc- In the lower figure, the solid lines represent the solar density p{x) and 
a linear approximation to p{x) that has the correct initial and final densities and the 
correct slope at Xc- (Thus the linear and exact density would almost exactly correspond 
over the nonadiabatic region illustrated in the upper figure) . The MSW equation can be 
solved analytically for the hnear wedge. By extending the wedge to ±00 (dotted hues) 
and assuming adiabatic propagation in those regions of unphysical density, one obtains 
the simple Landau-Zener result of Eqs. (31) and (32). 

Figure 14: MSW conversion for a neutrino produced at the sun's center. The upper 
shaded region indicates those Srn^ /E where the vacuum mass splitting is too great to be 
overcome by the solar density. Thus no level crossing occurs. The lower shaded region 
defines the — sin^2^^ region where the level crossing is nonadiabatic (7c < 1). The 
unshaded region corresponds to adiabatic level crossings and thus to strong — > 
conversion. 

Figure 15: The MSW solutions allowed at 95% c.l. by the combined results of the 
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Homestake, SAGE/GALLEX, and Kamiokande experiments, including Kamiokande II 
day-night constraints, given the BP flux predictions. From Hata (1995). 

Figure 16: MSW survival probabihties P^^'^{E) for typical small-angle {Srn^ ~ 6 • 
10-%V2,sin2 2^^ ~ 6 • IQ-^) and large-angle {Sm'^ - IQ-^eV^, sin^ 2^^ ~ 0.6) solutions. 

Figure 17: Schematic of the SNO detector now under construction in the Creighton 
#9 nickel mine, Sudbury. Provided by R.G.H. Robertson and J.F. Wilkerson (private 
communication) . 

Figure 18: Photograph of the Kamioka Mine cavity that will house Superkamiokande. 
Photo provided by the Institute for Cosmic Ray Research, University of Tokyo. 
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Table 1: Neutrino fluxes predict by the BahcaU/Pinsonneault (with and without He 
diffusion) and Turck-Chieze/ Lopez standard solar models. 

Source (MeV) flux (cm-^s-^) 

BP (with diffusion) BP (without) TCL 

p+p^^H + e+ + 1/ 
i3N^i3c + e+ + u 

i5o^i5N+e+ + 1/ 
i7F^i7o+e+ + ly 
8B^8Be+e+ + u 
^He+p^*^He+e+ + 1/ 
^Be+e~ ^^Li + u 



0.42 


6.00E10 


6.04E10 


6.03E10 


1.20 


4.92E8 


4.35E8 


3.83E8 


1.73 


4.26E8 


3.72E8 


3.18E8 


1.74 


5.39E6 


4.67E6 




~15 


5.69E6 


5.06E6 


4.43E6 


18.77 


1.23E3 


1.25E3 




0.86(90%) 


4.89E9 


4.61E9 


4.34E9 


0.38(10%) 








1.44 


1.43E8 


1.43E8 


1.39E8 
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Table 2: Predicted capture rates in SNU of the BP and TCL SSMs for the =^^C1 and 
SAGE/GALLEX experiments. 



capture rates 



neutrino source 


37C1(BP) 


37C1(TCL) 


7iGa(BP) 


7iGa(TCL) 


DD 


0.0 


0.0 


70.8 


71.1 


pep 


0.2 


0.22 


3.1 


2.99 


^Be 


1.2 


1.10 


35.8 


30.9 


«B 


6.2 


4.63 


13.8 


10.77 


13N 


0.1 


0.063 


3.0 


2.36 


150 


0.3 


0.21 


4.9 


3.66 


Total 


8.0 


6.36 


131.5 


122.5 
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